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LEGAL  NOTICE 


This  report  was  prepared  as  an  account  of  Government  sponsored 
work.  Neither  the  United  States,  nor  the  Department  of  Defense, 
nor  any  person  acting  on  behalf  of  the  Department  of  Defense: 

A.  Makes  any  warranty  or  representation,  express 

or  implied  with  respect  to  the  accuracy,  completeness 
or  usefulness  of  the  information  contained  in  this 
report,  or  that  the  use  of  any  information,  apparatus, 
method,  or  process  disclosed  in  this  report  may 
not  infringe  privately  owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to  the  use  of, 
or  the  damages  resulting  from  the  use  of  any 
information,  apparatus,  method,  or  process 
disclosed  in  this  report. 

As  used  in  the  above,  "person  acting  on  behalf  of  the  Department 
of  Defense"  includes  any  employee  or  contractor  of  the  Department, 
or  employee  of  such  contractor,  to  the  extent  that  such  employee 
or  contractor  of  the  Department  of  Defense,  or  employee  of  such 
contractor  prepares,  disseminates,  or  provides  access  to,  any 
information  pursuant  to  this  employment  or  contract  with  the 
Department  of  Defense,  or  his  employment  with  such  contractor. 
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I.  INTRODUCTION 

This  is  Volume  I  of  the  second  semi-annual  and  final  technical 
summary  report  describing  progress  of  the  research  and 
development  program  on  thermionic  conversion  of  heat  to 
electricity  under  the  Department  of  the  Navy,  Bureau  of  Ships 
Contract  NObs -90496.  The  work  performed  under  this  contract 
is  based  upon  and  is  a  continuation  of  work  initiated  under 
contract  NObs -88578.  Both  contracts  were  issued  by  the  U.  S. 
Navy,  Bureau  of  Ships,  and  supported  by  the  Advanced 
Research  Projects  Agency. 

The  objective  of  this  program  is  the  development  of  a  nuclear 
thermionic  electrical  generating  system  for  the  Department  of 
the  Navy.  Work  under  this  contract  is  being  performed  to 
develop  the  material  capabilities  which  are  essential  for  this 
nuclear  thermionic  system.  This  work  consists  of  three  major 
tasks.  They  are: 

TASK  A  -  UQ,  -  Refractory  Metal  Reaction  Kinetics 
This  task  includes: 

(1)  A  continuation  of  a  series  of  1000  hour  diffusion  ex¬ 
periments  of  U02-molybdenum  and  UO2 -tungsten 
fuel-emitter  pairs  in  the  temperature  regime  of 
1600  to  1800°C  for  molybdenum  and  1800  to  2200°C 
for  tungsten.  The  purpose  of  this  effort  is  the  quan¬ 
titative  identification  of  uranium  diffusion  (if  any) 
through  the  refractory  metal  and  the  limits  of  the 
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TASK  B 


TASK  C 


uranium  diffusion  rates.  This  will  be  accomplished 
by  alpha  particle  counting  techniques.  Metallo- 
graphic  and  electron  microprobe  analysis  will  be 
employed  to  document  the  extent  and  nature  of  solid 
state  reactions. 

(2)  The  determination  of  the  diffusion  rate  of  transport 
of  oxygen  through  molybdenum  or  tungsten  in  the 
temperature  range  of  1500  to  1800°C.  Thin  walled 
tubes  of  refractory  metals  will  be  heated  electrically 
in  a  cover  gas  of  flowing  purified  argon.  Known 
particle  pressures  of  oxygen  will  be  introduced  in¬ 
side  the  metal  tubes  and  the  argon  steam  monitored 
for  oxygen  content. 

-  Insulator  Materials  Development 

This  task  includes  the  identification  of  the  gaseous  species 
and  products  which  exist  when  high  purity  alumina  bodies 
and  niobium  and/or  niobium- 1%  zirconium  are  heated  to 
temperatures  of  interest  to  elucidate  possible  gas -solid 
reactions  in  a  thermionic  converter. 

-  Ceramic-to-Metal  Seal  Development 
This  task  includes: 

(1)  The  development  of  cesium-resistant  metallizing 
coatings  for  pure  alumina  ceramics. 
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(2)  The  development  of  molybdenum-alumina  cermet 
structures  for  use  as  emitter -collector  seal  in¬ 
sulators. 

(3)  The  evaluation  and/or  development  of  methods  of 
bonding  the  metallic  surface  of  insulators  to  struc¬ 
tural  metal  members. 

(4)  Long-time  testing  of  materials  and  seals  at  high 
temperatures  in  both  vacuum  and  cesium  vapor 
environments. 


The  work  accomplished  under  Tasks  A  and  B  is  included  in  this 
volume.  Task  C  is  reported  in  Volume  II. 


The  over -all  project  responsibility  for  contract  NObs -90496  at  the 
General  Electric  Company  is  held  by  Dr.  John  E.  VanHoomissen, 
Manager,  Nuclear  Thermionic  Power  Operation.  The  technical 
contributors  with  primary  responsibility  for  the  individual  tasks 
are: 


TASK  A-l 
TASK  A -2 

TASK  B 

TASK  C 


A.  I.  Kaznoff,  Nucleonics  Laboratory 

R.  A.  Ekvall,  A.  I.  Kaznoff,  L.  N.  Grossman, 
Nucleonics  Laboratory 

L.  N.  Grossman  and  A.  I.  Kaznoff,  Nucleonics 
Laboratory 

R.  H.  Bristow,  Tube  Department 
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II.  SUMMARY 

The  work  performed  under  contract  NObs -90496  is  described  in 
detail  in  Sections  III  of  this  volume  and  Volume  II.  Volume  II 
also  presents  in  Section  II  a  summary  of  the  work  accomplished  in 
Task  C.  Tasks  A  and  B  are  summarized  below. 

TASK  A  -  UC>2  -  Refractory  Metal  Reaction  Kinetics 
Sub -task  A-l  -  UO^  Diffusion  Studies 

The  diffusion  rates  of  uranium  through  tungsten  when  the 

tungsten  is  used  to  encapsulate  urania  was  determined  in 

the  temperature  range  between  1800  and  2000°C.  Similarly 

the  diffusion  rates  of  uranium  through  molybdenum  when  the 

molybdenum  is  used  to  encapsulate  urania  was  determined 

in  the  temperature  range  between  1600  and  1750°C.  The 

tests  were  performed  with  vapor  deposited  tungsten  and 

slightly  hyperstoichiometric  urania  and  with  several  types 

of  molybdenum  using  both  hyper -and  hypostoichiometric 

urania.  The  measured  diffusion  rates  with  the  urania- 

_  2  o 

tungsten  were  6.  33  x  10  gm/hrcm  at  1800  C  and 
<2  x  10  ^  gm/hrcm^  at  2000°C.  With  urania-molyb- 
denum  the  rates  ranged  from  1.  1  x  10  ^  gm/hrcm^  at 
1600°C  to  4  x  10  ^  gm/hr  cm^  at  1750°C  and  the  rates 
varied  with  the  type  of  molybdenum  and  the  stoichiometry 
of  the  urania.  It  was  generally  concluded  that  the  steady 
state  diffusion  rates  arc  of  such  magnitude  that  uranium 
diffusion  through  both  tungsten  and  molybdenum  is  not 
likely  to  affect  their  performance  as  thermionic  emitters 
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Sub-tas k  A  -2  -  Oxygen  Diffusion  Studies 

The  activation  energy  and  diffusivities  for  the  diffusion 
of  oxygen  through  polycrystalline  molybdenum  were 
determined  in  the  temperature  range  1500  1o  1800°C. 

The  two  experimental  technique's  used  to  measure  the 
diffusivities  yielded  comparable  results.  '1  he  activation 
energy  and  diffusivities  were  based,  in  part,  on  inter¬ 
polation  of  existing  data  on  the  solubility  limit  of  oxygen 

in  molybdenum.  The  activation  of  energy  was  67.2  +  20°'i 

-  8 

Kcal/mole  and  the  diffusivities  ranged  from  5  x  10 

l  i— *  n 

cm  /sec  at  1500°C  to  7.8x10  1  cm  /sec  at  1800°C.  The 
diffusivities  are  known  to  within  +50%. 

TASK  B  -  Insulator  Materials  Development 

A  mass  spectrometric  and  thermogravi  met  ric  study  of 
the  equilibrium  interaction  between  alumina  and  niobium 
between  1500  and  1900°C  has  been  conducted.  The  reaction 
proceeds  principally  as  A^O-^  +  3Nb  2Al(g)  +  3Nb0(s); 
a  secondary  reaction  of  almost  equal  importance  under 
neutral  conditions  is  Al^O-j  +  Nb  +  1/2  0^  A 1  ^ 0 ( fg )  + 

NbC^Cg)  +  Nb0(s).  The  latter  reaction  can  dominate  under 
oxidizing  conditions.  The  weight  loss  due  to  the  gaseous 
reaction  products  is  about  one  hundred  times  the  weight 
loss  due  to  evaportaion  of  alumina  under  neutral  conditions. 
Calculations  show  that  some  heavier  stable  oxides  (e.g., 
yttrium  oxide  and  thorium  oxide  will  not  interact  signifi- 
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cantly  with  niobium.  Gaseous  reaction  products  over 
beryllia -niobium  are  comparable  to  the  decomposition 
pressure  of  beryllia. 
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m.  TECHNICAL  PROGRESS 

TASK  A  -  UO-,  -  Refractory  Metal  Reaction  Kinetics 
Sub-Task  A  - 1  -  UO.,  Diffusion  Studies 

Introduction 

The  incorporation  of  thermionic  converters  in  a  nuclear 
reactor  core  requires  the  use  of  emitter  structures  contain¬ 
ing  nuclear  fuel.  The  fuel -cladding  composite  structure 
operates  at  high  temperatures  (>1  500°C).  Selection  and 
specification  of  materials  to  meet  these  needs  requires  the 
prior  determination  of  compatibility  of  the  fuel  with  the  clad¬ 
ding,  the  transport  rates  of  fuel  material  through  the  clad¬ 
ding,  the  contamination  of  the  cladding  surface  by  the  fuel 
components,  and  the  effects  of  thermal  cycling.  This  work 
was  directed  to  the  study  of  these  phenomena,  under  well- 
defined  conditions,  for  the  urania-tungsten  and  urania-molyb - 
denum  systems.  These  systems  represent  some  of  the  most 
stable  and  compatible  combinations  of  nuclear  fuel  and  emitter 
materials.  During  the  last  five  years  the  compatibility  of 

these  systems  was  established  by  exhaustive  experimental 
1  2 

work.  '  The  stability  of  urania  with  tungsten  and  molyb¬ 
denum  was  shown  for  periods  up  to  10,000  hours  at  thermionic 
design  temperatures. 

The  specific  objective  of  Sub  Task  A-l  was  to  continue  prior 
studies  on  the  diffusion  of  uranium,  from  urania,  through 
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tungsten  and  to  initiate  studies  for  molybdenum.  The  tests 
for  tungsten  were  to  be  performed  in  the  range  of  1800  - 
2200°C  and  for  molybdenum  primarily  in  the  range  of  1600  - 
1700°C.  The  magnitudes  of  diffusion  fluxes  were  to  be 
ascertained  by  alpha-counting  techniques;  auxiliary  analytical 
techniques  were  to  be  employed  to  verify  the  findings. 

Expe r imental  Techniques  &r  Procedures 

The  study  of  interactions  of  fuel  with  the  cladding  material 
was  conducted  with  encapsulated  fuel  (Figure  1).  The  dif- 
fusional  transport  was  observed  across  the  20-mil  thick 
"window"  in  the  capsule.  The  provision  of  much  thicker 
walls  and  end -closure  served  to  limit  the  diffusion  transport 
in  directions  other  than  through  the  window.  Consequently, 
the  capsule  system  approximated  a  one -dimensional  diffusion 
experiment.  The  capsule  was  heated  by  electron  bombard¬ 
ment  on  the  capsule  face  opposite  the  window.  This  method 
of  heating  produced  temperature  gradients  in  the  capsule  and 
in  the  fuel.  The  temperatures  were  measured  by  optical 
pyrometry  on  the  two  black-body  holes  in  the  thick  walls 
of  the  capsule.  This  pyrometry  established  the  temperature 
at  the  cladding  "window"  -  fuel  interface  and  the  gross 
temperature  gradient.  The  temperature  gradient  across  the 
fuel  was  calculated  to  be  200°C/inch  at  1800°C  and  480°C/inch 
at  2200°C,  The  test  temperatures  cited  refer  always  to  the 
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Overall  Height 
Sate  to  Shoulder 
Diameter 

Window  Thickness 


0.340  inch 
0.265  inch 
0.435  inch 
0.020  inch 


FIGURE  1.  DIFFUSION  CAPSULE 
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Assembly 


9 
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fuel-cladding  interface  at  the  diffusion  "window".  The 
presence  of  temperature  gradients  simulated,  qualitatively, 
the  thermal  regime  in  the  reactor  where  temperature  gradi¬ 
ents  are  always  present.  These  temperature  gradients  are 
of  major  importance  because  they  provide  the  driving  force 
for  urania  redistribution,  i.  e.  ,  the  sublimation  of  the  nuclear 
fuel  and  its  redeposition  at  regions  of  lower  temperatures 
(the  cladding).  Deposition  on  the  cladding  provides  and 
ensures  intimate  contact  between  fuel  and  cladding  and  sub¬ 
stantially  lowers  the  thermal  impedance  of  the  fuel-cladding 
interface.  Temperature  gradients  can  also  account  for  the 
sublimation  fractionation  of  the  more  volatile  phases  formed 
in  the  fuel  body  as  a  result  of  stoichiometry  changes  or  fission 
product  buildup. 

-  7  -  8 

The  tests  were  made  in  vacuum  between  3x10  and  8x10 
torr.  The  diffusing  species  through  the  20-mil  cladding  thick¬ 
ness  were  free  to  volatilize  from  the  outer  surface.  The  con¬ 
densibles  of  the  diffusing  flux  were  trapped  on  cold  (20°C) 
aluminum  targets  which  were  placed  directly  opposite  the 
"window"  area  in  a  water-cooled  magazine.  The  basic  com¬ 
ponents  of  the  test  device  are  shown  in  Figure  2.  Periodically, 
the  targets  were  removed  through  a  vacuum  interlock  without 
interrupting  the  test  and  analyzed  for  uranium  content.  The 
capability  of  sampling  without  interrupting  the  test  conditions 
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was  a  very  desirable  feature.  Apart  from  increasing  the 
utilization  of  the  device  by  this  sampling  technique,  it  was 
also  noted  that  shutdowns  (cooldowns  and  heatups)  could 
affect  the  results. 

The  fuel  used  in  this  work  was  fully  enriched  urania.  This 
selection  was  made  as  a  result  of  two  considerations: 

(a)  To  differentiate  the  uranium  obtained  as  a  result  of 
diffusion  from  the  trace  residual  natural  uranium  found  in 
refractory  metals  and  (b)  To  increase  the  sensitivity  of  the 
alpha  counting  technique  by  virtue  of  the  high  U-2  34  and 
U-235  content  in  enriched  urania  which  is  not  available  in 
natural  urania. 

The  primary  analytical  technique,  alpha  counting,  showed 

that  high  sensitivity  could  be  achieved.  One  count  per  minute 

-8 

represented  1.  325  x  10  gm  of  total  uranium  transport. 

The  counting  was  sensitive  to  at  least  0.  1  count/min.  above 

.  _9 

background.  It  is  believed  that  a  sensitivity  to  1.  3  x  10  gm 
of  uranium  was  achieved.  The  main  sensitivity  was  obtained 
from  the  U-234  content. 

The  secondary  analytical  technique  was  dilution  mass  spec¬ 
trometry  in  which  an  internal  calibration  (spike)  of  U-233  was 
used.  Generally,  the  results  were  obtained  on  U-235  and 

U-238  content  (where  possible  U-234  was  also  determined). 
The  determination  of  isotopic  content  was  necessary  to 
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eliminate  doubts  that  the  uranium  activity  observed  was  due 
to  the  enriched  fuel  which  was  diffused  and  not  due  to  the 
residual  uranium  in  the  test  components.  It  was  found  that 
U-238:U-235  ratios  were  higher  than  those  in  the  enriched 
fuel.  This  indicated  the  "boil -out"  of  the  natural  uranium 
impurity  in  the  refractory  metals.  This  contribution  was 
far  too  small  to  affect  alpha  counting  techniques,  which  were 
dependent  primarily  on  U-234,  or  the  contents  of  U-235 
obtained  by  mass  spectrometry. 

The  correspondence  between  mass  spectrometry  and  alpha 
counting  was  good,  especially  for  uranium  contents  about 

_9 

10  gm.  The  typical  correspondence  obtained  on  actual 
test  targets  was: 


_ Amount  of  U-235  (gm) _ 

Test  Alpha  Counting  Mass  Spectrometry 

1  1.  89  x  10 '7  1.  75  x  10’7 

2  1.  3  x  10'9  0.  54  x  10'9 

Generally,  mass  spectrometry  gave  lower  results  than  alpha 
counting  which  suggests  that  the  geometric  factor  in  alpha 
counting  was  lower  than  anticipated,  or  more  probably,  that 
there  was  some  urahium  loss  in  the  chemical  separation 
steps  leading  to  mass  spectrometric  determinations. 
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Mass  spectrometric  analysis  of  blank  targets  showed  that 
U-235  content  was  less  than  10  ^  gm.  Every  target  was 
alpha-counted  before  use  to  eliminate  the  possibility  of 
accidental  contamination.  Repeated  checks  on  target  cross  - 
contamination  in  the  experimental  apparatus  showed  negative 
results.  A  comparison  of  alpha-counting  and  mass  spectro¬ 
metry  results  indicated  that  contamination  of  targets  by  sub¬ 
limates  of  tungsten  and  molybdenum  did  not  affect  the  deter¬ 
mination  of  uranium  by  alpha-counting  to  any  appreciable 
extent. 

Supplementary  details  on  the  experimental  apparatus  and 
analytical  procedures  have  been  reported.  ^ 

Materials  Characterization 

The  UO^  chemical  analysis  was  carried  out  on  both  powder 
and  sintered  product.  The  sintered  pellet  product  was 
approximately  95%  of  theoretical  density,  and  was  cut  to  size 
by  an  ultrasonic  cutter. 
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TABLE  1 

Characterization  of  Urania 


Chemical  Analysis 


UO2  Powder 
as -received  (ppm) 

UO2  Sintered 
Pellets  (ppm) 

Comments 

A 1 

<5.0 

25.  0 

Increase 

B 

0.  1 

0.  6 

Unchanged 

C 

100.  0 

5.  0 

Large  decrease 

Cd 

<0.2 

<  1.0 

Unchanged 

Co 

1.0 

- 

Cr 

6.  0 

20.  0 

Increase 

Ca 

0 

1.0 

Unchanged 

F 

<50.  0 

<  5.0 

Decrease 

Fe 

20.  0 

75. 0 

Pick-up 

Mg 

4.0 

10.  0 

Virtually  unchanged 

Mn 

2.  0 

4.0 

Virtually  unchanged 

Mo 

11.0 

10.  0 

Unchanged 

Na 

1.0 

- 

Ni 

10.  0 

50.  0 

Increase 

Pb 

<1.0 

<  1.0 

Si 

20.  0 

60.  0 

Increase 

Sn 

<1.0 

<  1.0 

V 

- 

<  15.  0 

O/U 

2.  04 

2.0042 

Significant  decrease 

Density 

10.25  -  10.5 

gm/cc 

(theoretical  10, 

,  97  gm/cc) 

Isotope  Analysis 

Batch  I 

Batch  II 

U-235 

93.  50 

% 

92.78  % 

U  -234 

1.065  % 

1.25  % 

U -238 

5. 385  % 

5.743  % 
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The  tests  were  conducted  on  the  following  cladding  materials: 
Tungsten 


W-SFVD 


Mo-PM 


San  Fernando  vapor  deposited  (hydrogen 
reduction  of  WF^). 

Molybdenum 

General  Electric  power  metallurgy  sintered 
stock. 


Mo-AC  General  Electric  arc  cast. 

Mo-1/2  Te  General  Electric  Mo-1/2%  Ti. 


The  typical  analyses  for  these  materials  are  shown  on  the 
following  page  in  Table  2. 
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TABLE  2 


Chemical  Analyses  of  Refractory  Claddings  (in  ppm) 


Element 

W-SFVD 

Mo -AC 

Mo-PM 

C 

22-24 

100-210 

60-160 

O 

6-60 

30 

<  20 

N 

6-20 

20 

20 

Hf 

<  100 

<100 

A1 

<10 

<  50 

<  50 

Cd 

<100 

<100 

Cr 

<10 

<  50 

<  50 

Cu 

<10 

<  50 

<  50 

Fe 

<10 

<  50 

<  50 

Ni 

<  50 

<  50 

Si 

<10 

<  50 

<  50 

Sn 

<10 

<  50 

<  50 

Zr 

<25 

<100 

<  100 

Cb 

<  100 

<  100 

<  100 

Ta 

<  250 

<200 

<200 

Zn 

<100 

<100 

Ca 

<  50 

<  50 

Cb 

<  100 

<  100 

Na 

<100 

<  100 

Ti 

<10 

Mo 

<10 

Major 

Major 
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Test 


1. 

2. 


The  largest  variations  between  lots  and  within  lots  are  in 
carbon  and  oxygen  content.  With  the  exception  of  W-SFVD, 
the  materials  were  obtained  in  rod  form,  which  was  machined 
into  the  capsule  configuration.  Vapor -deposited  tungsten  was 
obtained  in  sections  suitable  for  incorporation  as  the  capsule 
structure.  Closures  on  capsules  were  made  by  TIG  welding. 
Experimental  Results.  Urania -Tungsten.  The  summary  of 
the  tests  conducted  on  UC^  qq^  anc*  tungsten  (vapor  deposited 
capsules)  are  given  below. 


T  ungsten 
Type 

Temperature 

rc) 

Test 

Duration 

(hours) 

Average 
Uranium 
Flux  ^ 

(gm  ^hr-cm  ) 

No.  of 

Thermal 

Cycles 

W-SFVD 

1800 

1273 

6.  33  x  10"9 

1 

W-SFVD 

2000 

1565 

<  2  x  10"10 

1 

Urania -Molybdenum.  A  summary  of  test  results  is  given  in 
Table  3.  (molybdenum  test  number  is  marked  with  an  asterisk 
to  differentiate  from  the  tungsten  test  in  further  discussions). 
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TABLE  3.  Urania-Molybdenam  Test  Results 


Test 

No.  of 

Test 

Molybdenum 

Urania 

Temp 

Duration 

Average  Uraniurn  Flux 
(gm/hr  -  crrC) 

Thermal 

No. 

_ Type 

0:U  Ratio 

(°C) 

(hours ) 

Cycles 

1* 

AC  -  1 

2. 0042 

1750 

250 

X 

o 

1 

o 

1 

1450 

380 

Transient  with  high 
peak  rate 

2 

AC  -1 

2.  12 

1750 

48 

High  initial  rate  of 

4  x  10“6  decreasing 
rapidly  to  10*9 

3* 

AC-1 

1.91 

1750 

320 

s-6  x  10  ^ 

1 

1605 

Decreasing  to 

3  x  10"10 

4* 

AC-1 

2.0042 

1600 

1  100 

1.  1  x  10"10 

2 

1700 

400 

1. 6  x  10- 10 

5 

AC  -2 

2. 0042 

1750 

300 

4  x  10"6 

2 

1600 

700 

i 

o 

6 

PM 

2. 0042 

1600 

1500 

-  7  -  8 

10  decreasing  to  10 

3 

7* 

Mo-1/2  Ti 

2.0042 

1600 

1000 

1.2  x  10'10 

Shown  in  Figure  3  are  the  results  of  Test  4  and  Test  1  *.  The  results  of 
Tests  2  and  3  are  shown  in  Figures  4  and  5. 
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Discussion  of  Results 

The  results  cited  for  vapor -deposited  tungsten  (Test  No.  1) 

represents  the  average  flux  over  the  test  period.  The  actual 

-  8  2 

diffusion  flux  started  at  a  high  rate  of  8.  2  x  10  gm/hr  crn 

_Q 

and  leveled  out  during  the  last  700  hours  to  2.  75  x  10  gm/hr 
2 

cm  .  Within  the  limitations  of  this  experimental  program 
and  analytical  capability  the  results  show  that  GE  powder 
metallurgy  tungsten,  tested  in  the  previous  program  year, 
showed  the  lowest  uranium  diffusion  rate  followed  by  arc -cast 
tungsten  and  then  followed  by  the  vapor  deposited  tungsten. 
However,  it  should  be  noted  that  surface  contamination  due 
to  absorbed  uranium  atoms,  had  been  shown  to  be  only  about 
1/100  of  a  monolayer  for  tungsten  claddings  used.  It  should 
be  also  noted  that  in  all  instances  the  transport  rates  were 
very  low  and  remained  low  at  higher  temperatures,  which 
are  above  the  temperatures  necessary  for  systems  design 
and  application.  It  is  not  known  why  the  test  at  2000°C  pro¬ 
duced  lower  apparent  rates  than  that  at  the  lower  temperature 
of  1800°C.  Vapor-deposited  tungsten  showed  very  minor 
grain  growth  and  the  columnar  grain  structure,  oriented 
parallel  to  the  diffusion  path,  remained  virtually  unchanged. 

In  this  structure  there  was  essentially  one  grain  in  the  20-mil 
cladding.  This  microstructure  may  account  for  the  larger 
diffusion  rates  when  compared  to  the  other  types  of  tungsten. 
The  observed  trend  is  higher  diffusion  rates  for  larger  grain 
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sizes.  In  all  cases  the  grain  size  in  the  diffusion  "window  ' 
area  of  the  capsule  was  the  same  as  in  the  rest  of  the  capsule 
body  which  indicates  that  this  diffusion  cannot  be  expected 
to  change  the  microstructure  and  electron  emission. 

Cycling  seemed  to  promote  the  penetration  of  urania  into 

the  grain  boundaries  to  a  limited  depth  of  about  1  to  3  mils. 

Electron  microprobe  examination  was  used  to  identify  these 

uranium -bearing  deposits.  It  is  postulated  that  the  mechanism 

for  this  penetration  is  as  follows:  On  cooldown,  free  uranium 

metal  is  produced  by  the  reaction: 

2  UO„  -UO,  +  x  U. 

—s -  2-x  2  z-r — 

2 -x  2+x 

The  elemental  uranium  penetrates  the  grain  boundary.  On 
reheating  to  the  original  temperature,  the  uranium  in  the 
grain  boundary  can  reoxidize  by  the  diffusion  of  oxygen  from 
the  UO2  next  to  it.  The  reverse  reaction  takes  place 
x  U  +  UO,  -  2  UO, 

v  “  >)  2  "X, 

2  -X  2-x 

This  results  in  a  racheting  mechanism  which  can  lead  to 
failure.  This  attack  is  not  general  and  is  not  observed  at 
all  grain  boundaries.  This  leads  one  to  suspect  that  only 
those  grain  boundaries  which  happen  to  be  in  contact  with  a 
uranium  precipitate  are  affected.  The  kinetics  of  this  process 
are  undoubtedly  affected  by  cooling  and  heating  rates,  the 
diffusion  rates  of  uranium  in  the  cladding,  and  oxygen  dif - 
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fusion  in  uranium  and  urania.  At  temperatures  of  interest 
to  thermionics,  such  attack  is  very  minor.  This  attack  can 
be  expected  to  increase  in  severity  with  higher  temperatures 
and  in  environments  favoring  a  higher  extent  of  sub -stoichi¬ 
ometry  (i.  e. ,  reducing  conditions)  in  the  urania  fuel  and 
power  cycling. 

The  urania  fuel  was  hypostoichiometric  after  the  test  in  all 
test  capsules;  final  0:11  ratio  was  of  the  order  of  1.  99.  The 
hypostoichiometry  increased  with  increase  in  temperature  as 
is  evident  from  metallographical  examination.  Minor  amounts 
of  tungsten  precipitates  were  found  in  the  urania  fuel.  These 
precipitates  were  probably  formed  in  the  initial  heatup  period 
when  the  oxygen  content  of  the  fuel  was  still  high.  Pore 
migration  up  the  temperature  gradient  in  urania  was  observed 
in  all  specimens. 

The  results  of  tests  in  the  system  Mo-UC^  showed  that 
diffusion  rates  at  1600  to  1750°C  ^  were  comparable  to 
those  in  the  W-UC^  system  at  1800  to  2200°C  ^  .  The  high 
rates  in  Tests  5*  and  6*  may  be  attributed  to  the  high  carbon 
content  in  the  molybdenum  but  a  definitive  explanation  is  not 
available  at  this  time.  The  main  result  of  Test  1*  was  that 
it  demonstrated  that  cycling  again  produced  high  peak  fluxes 
as  was  shown  in  an  earlier  study  (See  Figure  3  W-UC^). 
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Test  2*  was  designed  to  explore  the  effects  of  excess  oxygen. 

It  was  clearly  shown  that  very  high  diffusion  fluxes  resulted 

and  decayed  rapidly  with  time  (Figure  4).  This  was  not 

unexpected  because  of  previous  studies  of  the  reactions 

7 

in  the  Mo  +  system.  The  main  feature  of  the 

reaction  is  the  formation  of  low  metling  uranium -molybdenum - 
oxygen  compound.  This  compound  attacks  and  penetrates 
grain  boundaries  of  molybdenum  cladding.  In  a  reducing 
environment  and  at  high  temperatures,  this  compound  is 
unstable  and  decomposes  by  losing  oxygen,  to  yield  finally 
urania  and  molybdenum.  The  evidence  for  this  reaction  can 
be  readily  seen  by  inclusions  of  urania  in  the  grain  boundaries 
of  molybdenum  and  molybdenum  in  the  urania  body.  The 
effect  on  the  diffusion  rate  of  uranium  through  molybdenum, 
under  these  excess  oxygen  conditions,  is  a  large  increase 
compared  with  the  results  for  stoichiometric  urania.  Loss 
of  oxygen  which  brings  the  system  closer  to  stoichiometry 
reduces  the  diffusion  rate.  Internal  pressures  developed 


in  this  capsule  were  sufficient  to  distort  the  structure. 


By  testing  an  0:U  ratio  of  1.  91  in  Test  3  ,  it  was  shown 
that  higher  diffusion  rates  will  be  observed  if  the  fuel  con¬ 
sisted  of  urania  and  uranium  at  the  test  temperature.  This 
test  confirmed  the  hypothesis  advanced  for  the  cycling 
phenomena  observed.  As  a  consequence  of  these  tests,  it 
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can  be  concluded  that  no  benefits  accure  from  the  use  of 
hypo-  or  hyperstoichiometric  urania  from  the  criterion  of 
lowering  the  diffusion  rates. 

At  present,  it  is  premature  to  conclude  that  arc -cast  molyb¬ 
denum  is  superior  to  powder  metallurgy  molybdenum.  It  is 
conceivable  that  control  of  interstitial  impurities  may  be  an 
important  factor,  particularly  if  one  considers  that  impurities 
such  as  carbon  can  produce  reducing  conditions  and  may  be 
areas  susceptible  to  attack  on  the  grain  boundaries.  It  was 
established  that  in  all  tests  the  fuel  became  hypostoichiometric 
as  in  the  case  of  UO^-W 

UOds)  -  UO-,  (s)  +  x°2  T 
^  <s-x  2 

The  oxygen  was  lost  either  as  elemental  oxygen  or  as  metal 
oxides,  by  diffusion  and  volatilization  from  the  capsule.  On 
cooldown  (typically  at  or  below  1600°C)  the  hypo-stoichiometric 
phase  is  unstable  with  respect  to  urania  and  uranium. 

2  UO,  -UO,  +  x  U 

x -  2  -x  2  x - 

2 -x  2 -x 

The  elemental  uranium  could  then  dissolve  in  the  cladding  and 
produce  higher  concentration  gradients.  These  gradients 
increased  with  time  to  give  the  high  peak  rates  of  diffusion. 
Depletion  of  the  uranium  precipitate  source  caused  the  decay 
of  the  diffusion  flux.  This  physical  picture  is  consistent 


26 


GEST  -2062 


with  the  recently  determined  phase  diagram  for  uranium  - 

4 

urania  at  high  temperatures.  It  should  be  further  noted 

that  temperature  gradients  should  move  the  more  volatile 
uranium  to  the  cladding  interface. 

The  effects  of  thermal  cycling  and  shut  down  are  again 
graphically  illustrated  on  a  run  on  UO^  encapsulated  in  arc- 
cast  molybdenum  (Figure  5).  This  test  clearly  indicates  that 
peak  increases  of  uranium  flux  can  be  expected  as  a  result 
of  lowering  of  temperature  or  a  shut  down  followed  by  test¬ 
ing  at  the  previously  established  temperature.  Figure  6 
illustrates  similar  effects,  but  in  this  case,  the  molybdenum 
is  a  powder  metallurgy  product.  A  common  feature  in  all 
the  tests  was  a  higher  initial  flux  at  the  beginning  of  the 
experiment  followed  by  a  decay.  This  appears  to  correlate 
with  expected  oxygen  loss  and  adjustment  of  fuel  stoichiometry. 

Some  of  the  microstructur al  changes  observed  arc  similar  to 

3  8 

those  discussed  for  tungsten  '  .  The  salient  differences 

between  the  Mo-UO-,  system  and  the  W-UO^  system  were 
that  grain  growth  of  molybdenum  in  the  diffusion  "window  ’ 
area  was  retarded  by  comparison  with  the  growth  observed 
in  the  body  of  the  capsule.  This  may  indicate  that  grain 
boundary  contamination  as  a  result  of  diffusion  retards  the 
grain  growth  in  that  area. 
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In  all  capsules,  the  urania  was  slightly  hypostoichiometric. 

The  capsule  of  Test  3*  showed  that  most  excess  uranium 
was  depleted  in  the  fuel  body.  This  means  that  uranium  can 
migrate,  dissolve  in  the  cladding,  and  diffuse  out  of  the  sys¬ 
tem  if  the  fuel  is  hypostoichiometric  and  is  two-phase  at  the 
test  temperature. 

Conclusions 

Several  conclusions  reached  are:  (1)  Low  diffusion  rates  of 
uranium  through  tungsten  can  be  expected  when  tungsten  is 
used  to  encapsulate  urania  between  1800  and  2200°C.  The 
same  observation  applies  to  molybdenum  between  1600  and 
1750°C.  (2)  Surface  contamination  of  tungsten  is  not  a 

serious  problem  since  the  coverage  amounts  to  approximately 
1/100  of  a  monolayer.  (3)  It  was  found  that  powder  metal¬ 
lurgy  tungsten  showed  the  lowest  diffusion  rate,  followed  in 
order  by  arc -cast  and  vapor -deposited  tungsten.  No  definitive 
choice  is  available  on  molybdenum  in  terms  of  choice  between 
arc -cast  and  powder  metallurgy  stock.  (4)  Thermal  cycling 
can  produce  higher  peak  diffusion  fluxes  of  uranium  in  molyb¬ 
denum  and  tungsten  capsules.  (5)  Hypo-  and  hyper  stoichiometry 
in  the  fuel  results  in  the  enhancement  of  diffusion  rates.  (6) 

It  is  not  likely  that  steady  state  fluxes  measured  will  affect 

g 

thermionic  emission  as  already  verified  in  UC^-W 
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Sub -task  A -2  -  Oxygen  Diffusion  Studies 
General 

The  rate  of  diffusion  of  oxygen  through  the  cladding  of 

a  molybdenum  clad  uranium  dioxide  thermionic  emitter 

module  may  be  important  in  determining  the  operating 

performance  of  the  thermionic  fuel  element  (TFE).  At 

temperatures  above  1600°C  experimental  evidence  has 
9-12 

shown  that  uranium  dioxide  becomes  substoichiometric 
yielding  oxygen  according  to  the  reaction 

UO^  ^^2  x  +  Y  (w^ere  x  >  °)*  (1) 

Slow  diffusion  of  this  oxygen  through  the  Mo  would  probably 
not  be  deleterious  to  the  operation  of  the  TFE.  Rapid 
diffusion  of  this  oxygen  through  the  cladding,  however, 
would  cause  the  formation  of  volatile  molybdenum  oxides  at 
the  emitting  surface.  The  oxides  could  interfere  with  TFE 
performance  by  remaining  as  a  gas  in  the  interelectrode 
space  or  by  depositing  as  a  solid  product.  The  degree  of 
substoichiometry  of  the  fuel  also  would  increase  by  the 
loss  of  oxygen.  On  subsequent  lowering  of  the  temperature, 
disproportionation  of  the  fuel  would  occur  according  to  the 
reaction 

U°2-x  ~TT  u°2  +  T  U  (where  x  >  o).  (2) 

The  free  uranium  produced  by  this  reaction  could  react  with 
the  molybdenum  cladding.  Because  of  the  importance  of  the 
rate  of  oxygen  diffusion  in  molybdenum,  the  present  study 
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was  undertaken.  The  diffusion  coefficients  (diffusivities ) 
of  oxygen  through  molybdenum  were  measured  in  the  anti¬ 
cipated  thermionic  reactor  operating  temperature  range, 
1500  to  1  800°C  . 

Two  experimental  techniques  were  used  for  measuring 
the  diffusivity,  D.  The  first  approach  involved  passing 
a  stream  of  oxygen  gas  over  a  molybdenum  rod  and 
analyzing  the  resulting  gradient  by  using  the  non-steady 
state  solution  to  Fick's  second  law  for  an  infinite  cylinder. 
Because  of  experimental  difficulties  with  this  method,  a 
second  approach  was  also  used.  In  this  method,  the 
steady  state  gradient  and  flux  of  oxygen  moving  through 
a  thin -walled -molybdenum  cylinder  were  measured.  The 
diffusivity  was  calculated  using  Fick's  first  law.  Final 
deter mination  of  diffusivity  values  required  the  estimation 
of  the  solubility  of  oxygen  in  molybdenum  as  a  function  of 
temperature  . 

Experimental  Procedure 

Two  methods  were  used  to  measure  D,  a  "saturation" 
method  and  a  "flux"  method.  Each  technique  is  discussed 
individually  below. 

"Saturation"  Method 

The  first  method  involved  the  determination  of  diffusivity 
by  measuring  the  concentration  gradient  of  oxygen  in  a 
cylindrical  molybdenum  rod  subjected  to  treatment  in 
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oxygen  and  comparing  this  gradient  with  that  predicted  by 

the  solution  of  Fick's  second  law  diffusion  equation  for  an 

13  14 

infinite  cylinder.  (Details  were  given  previously  ). 

Initial  plans  required  the  uniform  saturation  of  molybdenum 
rod  samples,  at  a  given  temperature,  to  an  oxygen  con¬ 
centration,  Cq,  in  equilibrium  with  the  maximum  permis¬ 
sible  external  pressure  o'  oxygen.  (Cq  is  also  the 
concentration  of  oxygen  at  the  surface  of  the  rod  at  any 
time  during  exposure  to  the  oxygen  atmosphere).  The 
saturation  treatment  was  to  be  followed  by  a  deoxidation 
treatment  in  hydrogen  for  a  prescribed  time.  The 
resulting  oxygen  concentration  gradient  was  then  to  be 
measured  in  the  partially  deoxidized  samples.  Calcula¬ 
tions  of  saturation  times  based  on  values  of  D  previously 

15,  16 

estimated  from  published  low  temperature  data 
indicated  that  this  procedure  could  turn  out  to  be  inordin¬ 
ately  time  consuming.  Because  of  the  possibility  of 
excessive  time  for  experimentation,  the  procedure  of 
saturation  followed  by  deoxidation  was  replaced  by  oxida¬ 
tion  treatments  only.  This  change  was  not  expected  to 
significantly  affect  the  anticipated  results. 

A  further  simplification  was  introduced  in  the  experimental 
procedure  at  this  point  which  made  it  possible  to  determine 
D  and  CQ  from  experiments  formerly  designed  to  determine 
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One  method  of  measuring  Cq  at  a  given  temperature  is 

to  heat  molybdenum  rod  samples  of  radius,  a,  in  the 

maximum  permissible  oxygen  pressure  for  different  lengths 

of  time,  t.  The  average  oxygen  concentration,  after 

a  given  time  of  exposure  can  then  be  measured  and  plotted 

as  a  function  of  time.  As  the  time  of  exposure  increases 

the  oxygen  concentration  should  asymptotically  approach 

a  maximum,  or  saturation  value ,  C  .  The  time  required 

o 

to  reach  C  is  called  the  saturation  time.  The  curve  of 
o 

C  vs  t  is  called  a  saturation  curve, 
av 


If  Cq,  CF,  (the  uniform  oxygen  concentration  in  the  samples 

2 

before  heating  in  an  oxygen  environment)  and  Dt/a  are 

known,  a  value  for  D  can  be  found  from  the  experimentally 

determined  curve  of  oxygen  concentration  versus  time  (i.e.  , 

the  saturation  curve)  be  taking  advantage  of  the  following 

theoretical  observation.  It  can  be  shown  that  each  curve  of 

the  theoretical  family  of  curves  obtained  by  plotting 

r  C(r)  -  C.  /  Cq  -  C.  versus  r/a  has  a  unique  value 

2 

of  Dt/a  and  the  area  under  each  curve  A  is  unique.  r  is 
the  radial  distance  from  the  center  of  the  rod;  C(r)  is  the 
oxygen  concentration  at  r  after  a  given  time,  t  .  Further, 


it  can  be  shown  that  a  given  area  is  related  to  C  of  a 

av 

particular  rod  by  the  following  expression  (see  appendix  A). 


This  relationship  permits  the  determination  of  D  using  a 
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method  of  successive  approximations.  This  method  is 
outlined  below. 

An  estimate  of  D  can  be  made  from  the  saturation  time 

indicated  by  the  experimental  saturation  curve.  The  value 

2 

of  Dt/a  can  then  be  calculated  for  each  experimental 
point  of  the  experimental  saturation  curve  and  the  area 

2 

under  the  theoretical  curve  corresponding  to  each  Dt/a 

value  computed.  Knowing  A  and  having  experimentally 

measured  C  and  C.,  a  theoretical  value  of  C  corres- 
o  i  av 

ponding  to  each  experimentally  determined  can  be 

found  from  equation  (3).  Plotting  the  theoretical  C&v 
values  vs  time  yields  a  saturation  curve.  This  theore¬ 
tical  curve  is  then  compared  with  the  experimental 
curve.  If  the  theoretical  curve  does  not  closely  approxi¬ 
mate  the  experimental  curve  a  different  value  of  D  is 
estimated  and  the  above  procedure  repeated  until  good 
agreement  is  obtained  between  the  theoretical  and  ex¬ 
perimental  saturation  curves.  The  value  of  D  associated 
with  the  best  fitting  theoretical  curve  is  the  D  for  oxygen 
in  molybdenum  at  that  temperature. 

To  determine  the  experimental  saturation  curve  at  1500 

and  1800°C ,  1/16  -inch  rods  of  molybdenum  were  heated 

in  a  mixed  argon-oxygen  atmosphere.  The  partial  pressure 

of  oxygen  at  a  given  temperature  was  kept  approximately 

an  order  of  magnitude  below  the  equilibrium  oxygen 

1  7 

pressure  for  the  formation  of  solid  MoO^  Oxygen  partial 
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pressures  used  were  7.  6  x  10~7  torr  and  7.  6  x  10-5 
torr  at  1500  and  1800°C,  respectively.  The  samples 
were  supported  vertically  within  a  water-cooled  copper 
concentrator  inside  a  vacuum  chamber  and  inductively 
heated  by  an  RF  coil  external  to  the  vacuum  chamber. 
Pressures  were  measured  with  a  McLeod  guage.  Flow 
of  the  gas  into  the  vacuum  chamber  was  controlled 
mechanically  with  a  micrometer  needle  valve.  For  a 
given  oxygen  pressure,  micrometer  settings  varied  only 
slightly  from  sample  to  sample.  Sample  temperature 
was  measured  with  an  optical  pyrometer  within  +  1%.  A 
black  body  hole  (0.  020-inch  diameter  x  0.  100-inch  in 
depth)  at  the  upper  end  of  the  sample  and  a  black  body  hole 
in  the  molybdenum  sample  support  were  viewed  for  temp¬ 
erature  measurements.  Temperature  uniformity  over 
the  length  of  the  sample  was  kept  within  20°C.  The  appar¬ 
atus  is  shown  in  Figure  7. 

Difficulties  in  maintaining  a  steady  temperature  at  1800°C 
eventually  required  the  use  of  1 /8-inch  diameter  samples 
at  that  temperature  to  provide  better  electrical  coupling 
between  samples  and  the  induction  heater. 

Rods  of  both  diameters  were  given  vacuum  (approximately 
10  torr)  anneals  prior  to  exposure  to  the  oxidizing  envi¬ 
ronment.  The  1/16  and  1/18-inch  diameter  rods  were 
heated  at  1900°C  for  7  and  28  hours,  respectively.  The 
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FIGURE  7. 


INDUCTION  HEATING  APPARATUS 


GEST -2062 


purpose  of  the  anneals  was  to  reduce  C.  to  a  low  value  and 
to  change  the  wrought  structures  of  the  rods  to  uniform 
recrystallized  structures.  The  post-anneal  grain  size  of 
1/16-inch  diameter  rods,  however,  was  found  to  be  quite 
variable.  Photomicrographs  of  the  samples  before  and 
after  annealing  are  shown  in  Figure  8.  The  chemical 
compositions  of  the  rods  after  annealing  are  presented  in 
Table  4  .  The  precisions  of  oxygen  and  carbon  analyses 
given  in  this  report  were  +  8  and  +  5%,  respectively.  Where 
more  than  one  sample  was  analyzed  for  a  given  element, 

an  average  value  and  the  individual  values  (in  brackets)  are 
given. 


^  *  V4^v  mciliuu 


While  continuing  the  above  method  for  measuring  the  diffusi- 
vity  of  oxygen  in  molybdenum,  a  second  technique  was  ori¬ 
ginated  when  scatter  was  found  in  the  data  to  be  used  for 
determining  the  experimental  saturation  curve.  According  to 
Tick's  first  law,  the  flux,  J  (gms /cm2 /sec ),  or  quantity 
per  second,  of  diffusion  matter  passing  normally  through 
an  area,  A(cm2),  under  the  action  of  a  concentration 
gradient,  dc /dx  (gms /cm4)  is  given  by 

'"Oxygen  analyses  of  samples  were  made  by  Luvak  Incorporated,  Newton 
ass.  Carbon  and  nitrogen  analyses  were  performed  a.  National  Research 
Corporation,  Cambridge,  Massachusetts. 
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TABLE  4. 


Chemical  Composition  of  Molybdenum  Materials 


Rods 


1  /16-inch  dia^ 


Element 


inch  diax  '  1  /8  -inch  dia' 

(After  Vacuum  Anneal) 

4.9  4.6 


_ /,  a  uue 

1/16  -inch  dia  ' 
_ (As  -received) 


Caps 


Analytical 

Technique 


0 

(3.7 

’,7. 7,4. 6, 3. 7) 

(6. 

0,  3.  8,  3. 

9)  (8.5, 

7.  6,  7. 

8) 

24 

22 

18 

10 

36 

c 

(22,  13) 

{11 

,11.7) 

(44, 

27) 

52 

211 

N 

< 

10 

< 

5 

< 

5 

n 

.d.(° 

) 

n.d. 

A1 

10 

30 

< 

10 

< 

1 

< 

1 

B 

< 

2 

< 

5 

< 

5 

< 

1 

< 

1 

Ca 

20 

< 

10 

< 

10 

< 

10 

< 

10 

Cd 

< 

20 

< 

5 

< 

5 

< 

10 

< 

10 

Co 

< 

2 

50 

< 

10 

< 

5 

< 

5 

Cr 

100 

100 

< 

10 

< 

10 

< 

10 

Cu 

10 

< 

10 

< 

10 

< 

1 

< 

1 

Fe 

50 

500 

< 

5 

10 

10 

Mg 

1 

5 

< 

0.  1 

< 

1 

< 

1 

Mn 

1 

1 

10 

< 

1 

< 

1 

Mo 

Major 

Major 

Major 

Major 

Major 

Na 

50 

< 

50 

< 

50 

< 

50 

< 

50 

Ni 

5 

< 

10 

< 

10 

< 

10 

<r. 

10 

Si 

100 

500 

< 

5 

10 

10 

Sn 

< 

1 

10 

10 

< 

5 

< 

5 

Ti 

20 

< 

10 

< 

10 

< 

10 

< 

10 

V 

10 

< 

10 

< 

10 

< 

10 

< 

10 

W 

< 

20 

< 

100 

<100 

<1 

00 

< 

100 

Zn 

< 

10 

< 

10 

< 

10 

< 

10 

< 

10 

Zr 

< 

2 

50 

50 

< 

5 

< 

5 

Vacuum 

fusion 

Conducto  - 

metric 

Kjeldahl 

Emission 
spcctrogra - 
,,  Phic 


a  -  Sintered  rod  from  General  Electric  Company 
b  -  Arc  cast  rod  from  Climax  Molybdenum  Company 
c  -  Sintered  molybdenum  from  Cleveland  Tungsten  Company 
d  -  Sintered  rod  -  source  unknown 

e  -  n.d.  =  no  determination  3Q 
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J  =  -Dj^-  . 
dx 


(4) 


Solving  (4)  for  D  for  a  one -dimensional  case  gives 

D  =  J. 

Ac 


(5) 


It  is  assumed  in  equation  (5)  that  diffusion  is  taking  place 
in  only  one  dimension  and  that  D  is  independent  of  compo¬ 
sition.  Both  assumptions  are  met  in  the  experimental 
method  described  below.  The  diffusion  flux  passes  through 
a  thin  walled  cylinder  (0.  020-inch  thick)  which  closely 
approximates  one  dimensional  flow  and  the  solubility  of 
oxygen  is  <  100  ppm  in  molybdenum  and  thererfore  D  is  inde 
pendent  of  composition. 


To  measure  D  at  a  given  temperature  using  the  above  rela¬ 
tionship,  thin  walled  sealed  molybdenum  capsules  containing 
a  fixed  partial  pressure  of  oxygen  were  heated  in  vacuum 
(approximately  10  ^  torr).  The  flux  was  found  by  measuring 
the  total  weight  loss  rate  per  unit  surface  area  of  the  capsule 
and  correcting  the  weight  loss  for  oxygen  loss  by  mass  spec - 
trometrically  determining  the  oxygen  content  of  evaporating 
species.  The  concentration  gradient  was  established  by 
measuring  the  oxygen  concentration  of  molybdenum  rods 
contained  inside  and  outside  the  capsule  after  the  system  had 
been  heated  for  a  lehgth  of  time  estimated  to  saturate  the  rods. 

Test  capsule  preparation  was  as  follows.  Several  1/16-inch 
diameter  molybdenum  rods  were  placed  in  ultrasonically 
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cored  holes  of  a  sintered  pellet  of  Mo  +  MoO^  (approximately 

equal  amounts  of  each)  0.55 -inch  long  by  approximately  0.4- 

inch  diameter.  The  purpose  of  the  pellet  was  to  provide 

a  partial  pressure  of  oxygen  within  the  capsule  equal  to  the 

equilibrium  pressure  for  the  formation  of  solid  molybdenum 
17 

dioxide.  Under  these  conditions  the  maximum  material 
solubility  of  oxygen  in  molybdenum  at  a  given  temperature 
could  be  attained  at  the  inner  capsule  surface  without  nucle- 
ation  of  a  second  phase.  The  pellet  was  then  sealed  in  a 
molybdenum  capsule  (0.45 -inch  od  x  approximately  0.6 -inch 
long  x  0.  020-inch  wall)  formed  by  inert  gas  welding  machined 
molybdenum  lids  to  a  molybdenum  tube  section.  |a  11 
molybdenum  parts  were  vacuum  fired  (approximately  10"^  torr 
at  1500°C  for  one  hour  prior  to  assemblyj.  A  leak  check 
of  the  capsules  with  a  helium  mass  spectrometer  leak  detector 
showed  the  capsules  to  be  leak  tight.  A  longitudinal  and 
transverse  section  of  a  typical  test  capsule  (after  400  hours 
at  1500°C)  is  shown  in  Figure  9.  Photomicrographs  of  the 
as -received  structures  of  the  tube  cap  and  rod  are  shown  in 
Figures  10  and  1 1 .  Chemical  analyses  of  the  components  are 
given  in  Table  4. 

Two  capsules  were  isothermally  tested  at  1500°C  and  two  at 
1800°C  to  measure  weight  losses  and  to  establish  the  equili¬ 
brium  concentration  of  oxygen  in  the  rods.  The  tests  were 
carried  out  in  high  temperature  resistance  heated  vacuum 
furnaces.  Temperature  was  regulated  within  +  1%  of  the 
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ETCHED 


2.5X 


a.  Transverse  Section 


ETCHED  2.5X 


b.  Longitudinal  Section 

IGURE  9.  TYPICAL  MOLYBDENUM  DIFFUSION  CAPSULE 

(NO.  3)  CONTAINING  PELLET  OF  MOLYBDENUM 
AND  MOLYBDENUM-DIOXIDE  AND  MOLYBDENUM 
RODS,  AFTER  400  HOURS  VACUUM  ANNEAL. 
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ETCHED  50X 

a.  Tube 


ETCHED  50X 

b.  Cap 


FIGURE  10.  AS-RECEIVED  STRUCTURES  OF  MOLYBDENUM  TUBE 
AND  CAP. 
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ETCHED 


50X 


FIGURE  11.  AS-RECEIVED  STRUCTURE  OF  1  / 16  INCH  DIAMETER  ARC 
CAST  MOLYBDENUM  RODS  USED  INSIDE  AND  OUTSIDE  OF 
DIFFUSION  CAPSULES 
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desired  test  by  tungsten /rhenium  thermocouple  control 
systems. 

Mass  spec trometric  studies  were  made  to  determine 
evaporating  species.  Two  additional  capsules  were  tested 
in  a  Bendix  Time -of-Flight  mass  spectrometer  and  the 
species  determined  at  each  temperature  of  interest. 

Results  and  Discussion 

"Saturation"  Method 

Results  of  tests  to  determine  experimental  saturation  curves 
at  1500  and  1800°C  are  given  in  Figures  12  and  1  3.  The  number 
adjacent  to  each  point  is  the  sample  number.  The  magnitude 
of  the  sample  number  indicates  the  relative  chronological 
occurrence  of  the  test,  e.g.  ,  sample  1  was  tested  first. 

Weight  loss  associated  with  each  sample  is  plotted  as  a  function 
of  time  in  Figures  14  and  15  for  1500  and  1800°C  tests,  respec¬ 
tively.  The  theoretical  weight  loss  of  molybdenum,  based  on 

Langmuir  evaporation,  i.e.,  the  maximum  rate  of  evaporation 

19 

into  a  vacuum  was  calculated  and  is  also  shown.  The 

_2 

calculated  loss  of  molybdenum  is  only  6x10  mg  after 
409  hours  at  15  00°C.  It  is  obvious  that  molybdenum  oxide  must 
also  have  been  forming  and  evaporating  during  tests  at  both 
temperatures.  This  mechanism  would  account  for  the  additional 
material  loss  observed.  The  most  likely  gaseous  specie  is 
molybdenum  dioxide.  Appreciable  reduction  in  sample  diameter 


Time-Hours 


FIGURE  12.  AVERAGE  OXYGEN  CONCENTRATION  OF  1/16  INCH 
DIAMETER  MOLYBDENUM  RODS  VS  TIME  AT  1500°C 
(Pop  =  7.6  x  10”7  torr) 
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Time-Hours 


FIGURE  13.  AVERAGE  OXYGEN  CONCENTRATION  OF  1/16  AND  1/8 
INCH  DIAMETER  MOLYBDENUM  RODS  VS  TIME  AT 
1800  C  (Po?  =  7.  6  x  IQ'5  torr) 
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did  not  occur  until  samples  had  been  heated  for  190  hours  or 
more  at  1500°C.  At  1800°C,  however,  the  diameter  of  sample 
28  tested  for  only  52  hours,  was  reduced  by  approximately 
0.010-inch  (approximately  17%  radial  dec rease ),  Figure  16. 

As  might  be  expected  from  the  structure  of  the  1/16-inch 

diameter  rods  after  their  vacuum  anneal  (see  experimental 

procedure),  the  1 /16-inch  rods  tested  at  1500  and  1800°C 

exhibited  highly  variable  grain  size  which  could  not  be 

correlated  with  trends  in  C  values  of  the  samples.  Etched 

av  r 

structures  of  rods  are  shown  in  Figure  17.  The  1  /8-inch 
diameter  rods  showed  huge  equiaxed  grains  after  annealing, 
some  covering  an  entire  transverse  section.  The  results 
obtained  at  each  temperature  are  discussed  below. 

1 500°C  Tests 

The  data  of  Figure  12  show  considerable  scatter,  but  the 
saturation  curve  determined  by  points  1,2,  3, 4,  and  23  was 
selected  for  analysis  to  evaluate  D.  This  curve  was  selected 
because  the  five  samples  named  were  tested  in  the  same 
experimental  unit.  Test  number  81,  carried  out  in  a  different 
apparatus,  also  falls  on  the  curve  based  on  the  five  most 
reliable  tests . 

Samples  1,2,3,4,23,  and  46  were  all  tested  in  one  vacuum 
induction  apparatus.  The  rest  of  the  samples  were  tested  on  a 
second  apparatus  different  from  the  first  only  in  one  respect; 
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the  second  unit  had  a  larger  pumping  speed  and  consequently 
a  higher  throughput  of  the  argon-oxygen  gas  mixture  to  main¬ 
tain  the  desired  partial  pressure  of  oxygen.  The  scatter  of 
all  data  obtained  after  the  initial  series  of  runs  (i.e.,  1,2, 
3,4,  and  23)  was  so  great  that  no  attempt  was  made  to  con¬ 
struct  a  curve  through  them.  The  catter  may  have  been  the 

result  of  an  acute  sensitivity  of  C  to  small  variations  in 

av 

pressure  and/or  flow  rate. 

The  method  of  successive  approximations  outlined  in  the 
experimental  procedure  was  applied  to  points  1,2,3,  and  4. 

The  C  (25  ppm)  of  sample  23  was  treated  as  C  and  the 
average  oxygen  concentration  of  rods  before  testing  (4.  9 
ppm)  was  used  as  C..  Calculated  curves  for  a  transfer 

coefficient,  p,  -  (see  appendix  A)  most  closely  approximated 

-  8 

the  experimental  results.  Of  these,  the  curve  for  D  4x10 
2 

cm  /sec  showed  the  best  fit.  Figure  18  displays  the  error 

band  (shaded  area)  of  the  analytical  oxygen  values  for  the 

experimental  points  and  curve  A  represents  the  calculated 

-  8 

curve  for  D  -  4  x  10  .  The  saturation  time,  Tg,  associated 

with  curve  A  was  approximately  43  hours.  Since  there  was  a 
possible  error  of  +_  8%  in  the  analytically  determined  value  of 
any  Cav,  Co  might  also  be  as  low  as  23  ppm.  The  calculated 
curve,  B,  using  the  value  of  Cq  yielded  a  better  fot  to  experi¬ 
mental  data.  Values  of  D  and  T  associated  with  curve  B  were 

s 

-8  2 

3  x  10  cm  /sec  and  53  hours,  respectively.  Using  the  ether 
experimental  extreme  value  of  CQ  =  27  ppm  yielded  a  curve 
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which  did  not  fit  the  data  as  well  as  either  curve  A  or  B. 

One  can  expect  that  the  value  of  D  is  known  to  about  +  50%. 

1800°C 

Experimental  difficulties  at  1800°C  prevented  determination 

of  a  saturation  value  using  1/16 -inch  diameter  rods,  Figure 

13.  (Point  34  is  neglected  in  drawing  the  curve  for  1/16- 

inch  diameter  samples  because  comparison  with  1500°C  data 

shows  it  to  be  unrealistically  low).  Results  with  1  /8-inch 

diameter  samples  were  no  more  promising  and  no  value  of  D 

could  be  measured.  The  C  values  for  the  1  /8-inch  rods 

av 

were  quite  low  even  though  longer  saturation  times  were 

required  for  these  samples.  The  weight  loss  data  for  the 

1 /8-inch  rods  also  seem  low  relative  to  the  l/16-inch 

diameter  rods.  These  low  values  of  C  and  weight  are  not 

av  ° 

as  yet  understood. 

"Flux"  Method 

Limited  success  in  determining  D  using  experimental  method  I, 
led  to  the  origination  of  method  II  described  in  the  experimental 
procedure.  This  method  involved  the  experimental  evaluation 
of  the  various  parameters  related  by  Fick's  first  law.  Results 
of  vacuum  anneals  at  1500  and  1800°C  are  shown  in  Table  5. 

The  annealing  times  were  based  on  a  very  conservative  estimate 
of  the  time  it  would  take  at  a  given  temperature  to  saturate  the 
1/16 -inch  molybdenum  rods  inside  the  capsules.  The  weight 
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loss  rate  was  measured  by  weighing  each  capsule  after  each 
of  a  series  of  isothermal  short  term  (1/2  to  2  hours)  rund 
during  the  early  stages  of  the  vacuum  anneal.  The  average 
weight  losses  calculated  from  these  measurements  are  given 
in  Table  5.  If,  following  the  short  term  runs,  the  weight  loss 
rate  for  each  capsule  was  maintained  during  the  remainder 
of  the  isothermal  anneals  the  expected  weight  loss  would  have 
been  that  shown  in  Table  5.  The  actual  weight  losses  for  the 
1500°C  capsules  were  close  to  the  expected  losses  and  indi¬ 
cated  that  oxygen  was  available  inside  the  capsule  during  the 
entire  vacuum  anneal.  The  low  actual  weight  losses  for  the 
1800°C  capsules  indicated  a  loss  of  all  the  oxygen  long  before 
the  completion  of  the  1800°C  annealing  treatments.  X-Ray 
diffraction  and  metallographic  analysis  confirmed  the 
expectations  based  on  the  weight  loss  data.  Molybdenum  dioxide 
was  found  in  both  1500°C  capsules  but  only  molybdenum  was 
found  in  the  1800°C  capsules  after  testing.  In  the  1800°C 
capsule  only  a  sponge  of  molybdenum  was  left  where  once 
there  had  been  a  pellet  containing  MoC^  +  Mo.  Macrophoto¬ 
graphs  of  capsules  No.  3  and  1  are  shown  in  Figures  9  and  19. 
respectively.  The  variable  grain  size  of  the  molybdenum  rods 
tested  at  1500°C  and  the  very  large  grains  of  the  molybdenum 
rods  at  1800°C  are  evident.  Portions  of  the  caps  and  walls 
of  each  capsule  are  presented  in  Figure  20  and  21.  It  was 
assumed  in  the  concentration  gradient  calculation  that  no  sig¬ 
nificant  change  occurred  in  the  wall  thickness  of  the  capsules 
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FIGURE  20.  LONGITUDINAL  SECTIONS  OF  CAPSULE  NO  3  AFTER 
409  HOURS  AT  1500°C. 
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FIGURE  21.  LONGITUDINAL  SECTIONS  OF  CAPSULE  NO.  1  AFTER 
54  HOURS  AT  1800°C. 
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during  the  initial  test  period,  during  which  time  the  weight 
losses  were  being  determined. 

The  oxygen  content  of  the  rods  inside  capsule  No.  2  is  obviously 
below  the  saturation  value  since  all  oxygen  escaped  from  the 
capsule  long  before  the  end  of  the  test.  The  oxygen  content  of 
the  rods  outside  of  capsules  No.  3  and  4  are  not  accurate 
since  a  coating  containing  tantalum  (high  solubility  for  oxygen) 
was  found  on  parts  of  the  rods.  The  tantalum  apparently  eva¬ 
porated  from  a  surrounding  support. 

Two  sealed  capsules  nearly  identical  to  those  used  in  the 
vacuum  anneal  studies  were  examined  mass  spectrometrically 
to  discover  the  emanating  vapor  species  between  1500  and 
1800°C.  Molybdenum,  molybdenum  oxide,  molybdenum  dioxide, 
and  molybdenum  trioxide  were  present  up  to  1700°C.  Data 
above  this  temperature  was  not  obtained  because  the  two 
capsules  ruptured  during  the  tests.  The  average  composition 
of  the  escaping  oxygen  species  was  MoOg  ^  +  q  5'  ^he  dominant 
specie  was  molybdenum  dioxide. 

The  flux,  J,  was  computed  for  capsules  1,  2,  3,  and  4  from 
the  known  weight  loss  rate,  evaporating  species  and  the 
external  geometric  capsule  surface  area. 

The  concentration  gradient  for  the  i500°C  capsules  was  calcu¬ 
lated  assuming  that  the  rods  outside  the  capsules  would  have 
reached  the  same  oxygen  concentration  as  those  outside  the 
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1800°C  capsules  if  tantalum  had  not  masked  their  true 

values.  Because  of  the  loss  of  oxygen  from  the  1800°C 

capsules,  a  value  for  the  oxygen  concentration  at  the 

inside  surface  of  the  capsule  wall  for  this  temperature 

was  estimated  to  allow  calculation  of  the  concentration 

gradient.  (Appendix  B).  Diffusivity  values  calculated 

from  Fick's  first  law  using  the  flux  and  concentration 

-  8 

gradient  data  discussed  above  were  5.  34  x  10  and 
4.  79  x  10  **  cm^  /sec  at  15  00°C  and  8.  06  x  10~^  and 
7.46  x  10  '  cm  /sec  at  1800  C.  The  error  in  the 
diffusivity  values  may  be  as  large  as  +_  50%. 

The  activation  energy  for  diffusion  was  determined  to 
be  67,200  cal/mole  from  the  slope  of  In  D  vs  1/T,  Figure 
22,  The  activation  energy  is  known  within  +_20%.  Values 
of  D  at  1600  and  1700°C  were  determined  from  Figure  16  to 
be  1.4  x  10  and  3.  5  x  10  cm  /sec,  respectively.  The 
diffusivity  appears  to  increase  by  about  a  factor  of  2  with 
every  100°C  increase  of  temperature.  This  change  is 
relatively  small  and  indicates  that  it  would  be  difficult  to 
obtain  reliable  data  at  temperature  intervals  <100°C. 

To  obtain  supplementary  values  of  D  in  the  1400  to  1800°C 

11 

range,  evaluation  of  data  obtained  by  Grossman,  et.al. 
was  extended.  In  determining  the  compatibility  of  a  capsule 
system  of  molybdenum  clad  UO2  135*  Grossman  measured 
weight  loss  as  a  function  of  temperature  as  shown  in  Table 
6.  The  weight  loss  was  caused  by  the  escape  of  excess 
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of  oxygen  from  the  capsule.  Grossman  found  thai  temperatures 
in  excess  of  1700°C  were  necessary  to  cause  significant 
interaction  between  UO£+x  and  molybdenum.  Only  trace 
amounts  of  molybdenum  were  found  along  the  periphery  of  the 
U02+x  pellet  of  a  capsule  tested  to  1700°C.  The  concentra¬ 
tion  of  oxygen  at  the  inside  capsule  surface  was  believed  to 
be  close  to  the  equilibrium  value  becaus  molybdenum  oxide 
was  found  over  much  of  the  inside  surface.  Because  of  the 
similarity  of  Grossman's  study  to  the  present  work,  the 
weight  loss  data  of  one  of  his  capsules  tested  at  various  temp¬ 
eratures  up  to  1700°C  were  used  to  calculate  values  of  D  for 
oxygen  diffusion  in  molybdenum.  Assuming  the  average 
evaporating  specie  was  Mof^  ^  +  q  05’  wc^ht  l°ss  rates  and 
J  were  calculated  at  each  temperature.  Using  the  equilibrium 
values  of  the  oxygen  concentration  in  molybdenum  as  estimated 
in  Appendix  B,  the  diffusivity  of  oxygen  in  molybdenum  was 
calculated  for  each  temperature  and  is  shown  in  Table  III.  The 
value  for  D  at  1800°C  was  found  by  extrapolation  of  Grossman's 
data  plotted  as  In  D  vs  1/T  (Figure  22).  The  activation  energy, 
Q,  for  Grossman's  capsule  was  found  from  this  plot  to  be  59.6 
Kcal/mole. 


The  diffusivity  data  reported  above  are  summarized  and  com¬ 
pared  in  Table  7  with  other  diffusivity  results  on  the  molybdenum 

oxygen  system  found  in  the  literature.  In  studying  the  effect 

15 

of  impurities  on  the  welding  of  molybdenum,  Platte  used 
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20 

limited  data  from  a  Battelle  Memorial  Institute  report 

to  estimate  Q  at  1928°C  as  48,300  cal/mole.  In  his 

calculation,  Platte  assumed  the  frequency  factor,  D^, 

2 

equaled  1  cm  /sec.  In  a  strain  aging  study  on  molybdenum, 

16  21 

Brock  used  the  Cottrell  equation  based  on  the  Portevin- 

-2  2 

LeChatelier  effect  (assuming  Dq  -  10  cm  /sec)  to  estimate 
a  Q  of  53,  500  cal/mole.  The  diffusivities  shown  in  Table  IV 
for  each  author  were  calculated  by  solving  the  arrhenius 
relationship  for  D  assuming  Q  was  temperature  independent 
and  using  each  author's  selected  value  of  Dq. 

Stringer  and  Rudman  have  been  studying  the  diffusivity  and 

22 

solubility  of  interstitials  in  molybdenum.  Their  early  results 

show  a  diffusivity  of  10  ^  at  1700°C  and  an  activation  energy 

_2  2 

of  54,500  cal/mole  (assuming  Dq  =  10  ‘  cm  /sec)  for  oxygen. 
Because  of  the  preliminary  nature  of  this  data,  they  concluded 
their  agreement  with  Brock's  values  at  1700°C;  "can  only  be 
regarded  as  an  encouraging  coincidence.  " 

The  diffusivities  approximated  from  the  literature  are  shown 
to  indicate  the  order  of  magnitude  of  D  one  might  expect  in  the 
temperature  range  of  interest  and  to  demonstrate  the  degree  of 
uncertainty  in  the  previous  knowledge  of  D.  Neither  Platte  nor 
Brock's  studies  were  primarily  addressed  to  the  diffusion  of 
oxygen  in  molybdenum  and  they  only  infer  and/or  postulate 
the  diffusion  parameters.  The  work  of  Stringer  and  Rudman  is 
specifically  aimed  at  measuring  the  diffusion  parameters,  but 
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conclusions  based  on  their  early  results  would  be  premature. 

A  point  adding  some  confusion  to  the  results  of  the  above  authors 

is  the  value  of  Dq  chosen  for  calculating  Q.  For  example, 

-2  2 

had  Platte  used  D  =  10  cm  hse c  to  calculate  Q,  he  would 
o 

have  obtained  an  activation  energy  of  28,  150  cal/mole.  Powers 

23 

-2  2 

and  Doyle  measured  D  as  10  cm  /sec  in  the  group-V 

24  ° 

metals,  while  Jacobs  found  a  value  of  1.3  for  for  oxygen 
diffusion  in  tungsten.  Both  values  of  D  came  from  internal 
friction  measurements.  It  seems  likely  that  the  value  of  Dq 
for  group-VI  metals  would  be  greater  than  that  for  the  group  V. 

The  values  of  D  experimentally  determined  in  the  present  work 
by  the  flux  and  saturation  method,  fall  midway  between  those 
of  Platte  and  those  of  Brock  and  Stringer  and  Rudman.  The 
D  values  determined  from  the  present  work  are  considered  more 
reliable  than  those  estimated  from  published  literature  for 
a  number  of  reasons.  The  present  studies  were  specifically 
designed  to  measure  the  diffusivity  of  oxygen  through  molybdenum 
at  high  temperature.  The  experimentally  determined  values 
at  1500°C  are  in  very  good  agreement  even  though  they  resulted 
from  three  distinctly  different  types  of  experiments.  Also, 
good  agreement  was  obtained  at  1800°C  between  the  capsule 
experiments  of  the  present  work  and  Grossman's  data. 

Of  the  two  flux  experiments,  the  present  work  is  considered 
more  accurate  than  Grossman's  because  the  present  capsule 
tests  were  specifically  performed  to  measure  D  and  were 
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carried  out  under  conditions  most  closely  at 
a  steady  state. 

The  diffusivities  reported  above  were  actually  apparent 
bulk  diffusivities  in  that  the  contribution  of  grain  boundaries, 
if  any,  was  neglected.  The  rate  of  diffusion  along  grain 
boundaries  is  generally  more  rapid  than  through  the  bulk 
of  a  given  material.  Diffusivities  measured  in  polycry¬ 
stalline  material  thus  tend  to  be  higher  than  the  true  bulk 
value.  For  the  capsule  tests,  the  contribution  of  the  grain 
boundaries  to  the  diffusion  rate  was  probably  low,  if 
significant  at  all,  because  of  the  large  molybdenum  grains 
(ASTM  2-3)  and  the  high  test  temperatures  (i.e.,  at  high 
temperature  the  difference  between  rate  of  bulk  and  grain 

boundary  diffusion  decreases  and  the  bulk  contribution 

25 

tends  to  overpower  the  grain  boundary  contribution). 

The  accuracy  of  the  diffusivities  derived  from  the  capsule 
tests  is  somewhat  limited  by  the  small  number  of  samples 
tested.  As  mentioned  previously,  however,  the  calculated 
diffusivities  were  corroborated  by  the  results  of  at  least 
one  other  experimental  approach.  The  accuracy  of  the 
diffusivities  depend,  in  part,  on  the  correctness  of  the 
concentration  gradient  estimates  (Appendix  B).  More 
precise  values  of  the  rate  of  diffusion  of  oxygen  in  molybdenum 
may  be  possible  when  the  solubility  of  oxygen  in  molybdenum 
as  a  function  of  temperature  is  known.  Q  was  calculated  from 
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the  slope  of  the  curve  resulting  from  plotting  In  D  vs  1/T. 
The  accuracy  of  the  Q  value  thus  depends  on  the  accuracy 
of  D.  More  precise  values  of  D  should  therefore  yield 
better  values  of  Q. 

Conclusions 

The  diffusivities  and  activation  energy  for  the  diffusion  of 
oxygen  through  polycrystalline  molybdenum  were  determinec 
in  the  temperature  range  1500  to  1800°C.  These  values, 
which  depend  in  part  on  an  estimate  of  the  solid  solubility 
of  oxygen  in  molybdenum  as  a  function  of  temperature,  were 
67.  2  .+  20%  Kcal/mole  for  Q  and  the  diffusivities  were  as  sho- 
be  low: 

Temperature 


1500 

5.  1  x 

10"8  +  50% 

1600 

1.4  x 

10"7  + 50% 

1700 

3. 5  x 

10-7  +  50% 

1800 

7.8  x 

10"7  +  50% 

The  activation  energy  determined  is  somewhat  higher  than 
than  the  values  reported  in  the  literature. 

The  values  of  D,  however,  fall  roughly  midway  between  the 
values  estimated  from  dr  given  in  the  literature  and  are 
believed  to  be  more  reliable  than  the  literature  values. 


(cm^  /sec ' 
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TASK  B  -  Insulator  Materials  Development 
General 

For  a  metal  plus  oxide  system,  chemical  stability  can 
generally  be  assured  if  the  free  energy  of  formation  of 
the  oxide  is  more  negative  than  that  of  any  solid  oxide  of 
the  metal.  If  volatile  products  are  produced,  stability 
may  still  be  maintained  in  an  isothermal  closed  system. 
However,  in  an  open  system  or  in  a  temperature  gradient, 
a  reaction  yielding  a  volatile  product  is  favored  with  increasing 
temperature.  The  interaction  between  alumina  and  niobium 
is  a  reaction  which  yields  both  volatile  and  solid  products. 

Niobium  is  used  with  alumina  in  structural  and  in  ceramic  - 

to-metal  seal  applications  because  of  the  close  match  between 

their  thermal  expansivities.  Alumina-  niobium  seals  have 

been  developed  for  use  at  1200°C  for  nuclear  thermionic 

26 

converter  application.  The  long  term  stability  of  the 

alumina -niobium  combination  at  high  temperature  has  been 

in  doubt  due  to  the  chemical  similarity  between  niobium  and 

27 

tantalum.  Ackerman  and  Thorn  have  shown  that  tantalum 
enhances  the  evaporation  of  alumina  through  the  reaction 
AI2Q3  +  Ta  =  Al^Otg)  +  TaO^g)  plus  other  gaseous  products. 
They  point  out  that  at  high  temperatures,  the  weight  loss  of 
a  tantalum  cell  containing  alumina  was  essentially  the  same 
whether  or  not  an  effusion  hole  existed  in  the  cell;  it  can  be 
concluded  that  at  least  one  of  the  product  species  leaves  the 
cell  after  diffusing  through  the  tantalum  walls. 
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Experimental  Procedure 


28 

The  niobium-alumina  reaction  was  studied  via  the  Knudsen 
method.  Tungsten  effusion  cells  containing  an  equimolar 
mixture  of  niobium  and  alumina  were  electron  bombardment 
heated  as  shown  in  Figure  23.  The  Knudsen  cell  effluent  was 
directed  into  the  ionization  source  of  a  Bendix  Time-of- 
Flight  mass  spectrometer.  The  cell  chamber  was  valved 
from  the  spectrometer  to  facilitate  sample  change  without 
breaking  vacuum.  Temperature  capability  of  the  cell  heating 
system  has  been  demonstrated  to  2600°C. 


The  Knudson  cell  chamber  was  differentially  pumped  with 

two  8  1/s  ion  pumps  which  also  monitored  the  total  cell 

pressure.  The  pressure  in  the  source  region  of  the  mass 

spectrometer  was  monitored  with  a  Bayard -Alpert  type 

ionization  guage.  The  mass  spectrometer  vacuum  was 

maintained  by  a  mercury  diffusion  pump  which  was  baffled 

from  the  spectrometer  with  Freon-cooled  and  liquid  nitrigen- 

-8 

cooled  baffles.  Typical  source  pressure  was  2  to  8  x  10~ 

_7 

torr;  cell  chamber  total  pressure  was  between  2  x  10  and 
5  x  10'6  torr  during  experimentation. 


Sample  temperature  was  measured  by  sighting  into  the 
effusion  hole  with  an  optical  pyrometer.  The  pyrometer  was 
calibrated  against  an  NBS  secondary  standard  tungsten  strip 
lamp  between  900  and  2300°C.  Calibration  has  been  extended 
to  3000°C  by  the  rotating  sector  disk  technique.  Sample 
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FIGURE  23.  KNUDSEN  CELL  AND  ELECTRON  BOMBARDMENT  HEATER  ASSEMBLY 
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temperature  was  maintained  constant  to  within  +  5°C 
by  the  electron  bombardment  power  supply.  A  mock-up 
of  the  Knudsen  cell  and  its  heater  outside  the  mass 
spectrometer  revealed  no  detectable  temperature  gradi¬ 
ents  over  the  cell. 

Starting  materials  were  powders  which  were  pressed  into 
low  density  pellets  and  fired  at  1500°C  for  1  hour  at 
4  x  10  ^  torr.  The  starting  alumina  contained  less  than 
100  ppm  total  impurities  (principal  impurity  cations  were 
silicon  and  iron).  The  starting  niobium  contained  <400  ppm 
metallic  impurities,  but  had  the  following  interstitial 
impurity  content:  400  ppm  nitrogen,  1800  ppm  oxygen, 

27  ppm  hydrogen,  and  340  ppm  carbon.  The  oxygen 
content  of  the  starting  niobium  was  confirmed  at  4000  j_  3000 
ppm  by  X-Ray  determination  of  the  lattice  parameter;  the 
oxygen  was  presumed  to  be  in  solid  solution  in  the  starting 
metal  since  no  other  phases  were  detected. 

In  addition  to  the  elucidation  of  gaseous  products,  chemical, 
metallographic ,  and  X-Ray  analyses  were  made  of  the 
starting  materials  and  of  the  solid  residue  after  each  spec¬ 
trometer  or  weight  loss  experiment. 

Experiment  Results 

The  relative  abundance  of  various  mass  species  observed 
over  niobium  plus  alumina  is  shown  in  Figure24,  The 
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FIGURE  24.  MASS  PEAK  AMPLITUDE  VERSUS  TEMPERATURE  FOR  NEUTRALS 
ABOVE  Nb+Al,0, 
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common  gas  species  originated  in  the  cell  heater  assembly, 
not  inside  the  Knudsen  cell. 

The  mass  peak  intensities  shown  in  Figure  24  were  obtained 
after  3  days  continuous  heating  between  1200  and  1600°C. 

All  other  mass  peaks  were  at  most  one-tenth  as  intense  as 
those  shown  in  Figure  24  for  any  temperature.  No  significant 
changes  in  the  relative  intensities  were  observed  after  2  more 
days  of  heating  above  1500°C.  Both  neutral  and  ionized 
species  were  observed  entering  the  source  region.  Figure  25 
shows  the  temperature  dependence  of  the  ionized  species.  No 
other  ionized  species  were  observed.  The  relative  amounts 
of  ionized  and  neutral  species  changed  when  the  specimen 
position  was  changed;  the  ions  were  apparently  created  by 
collisions  between  neutrals  and  electrons  from  the  bombard¬ 
ment  filament.  No  corrections  to  the  observed  intensities 
have  been  made  for  ionization  efficiency  or  instrumental  fac¬ 
tors  . 

Residue  in  the  Knudsen  cells  was  analyzed  by  X-Ray  diffrac¬ 
tion  and  found  to  consist  of  three  phases.  Alumina  was  present 
with  no  change  in  its  lattice  parameter;  niobium  with  5000  ppm 

29 

oxygen  in  solution  (aQ  =  33  13  J?  4^  0.  005)  was  found;  and  a 

new  phase,  niobium  oxide,  was  present  (aQ  =  4. 2  12  j?  +_  0.  009). 

The  residue  was  examined  metallographically  and  the  three 

phases  were  seen.  To  bring  out  the  niobium  oxide  precipitates 

30 

in  niobium,  a  metallographic  procedure  after  Inouye  was 


followed. 
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FIGURE  25.  MASS  PEAK  AMPLITUDE  VERSUS  TEMPERATURE  FOR  IONS 
ABOVE  Nb+Al203 
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The  vapor  species  leaving  a  cell  containing  only  alumina 
were  also  examined  mass  spectrometrically.  The  species 
0(g),  Al(g),  and  02(g)  were  observed  in  the  ratios 
1:0.  5:0.  15  between  1570  and  2000°C .  No  corrections  for 
ionization  efficiency  or  instrumental  factors  have  been 
made  to  the  observed  ratios. 

The  observed  weight  loss  rate  of  a  tungsten  Knudsen  cell 

containing  niobium  plus  alumina  is  shown  in  Figure  26,  The 

data  in  Figure  24  inc hides  a  Clausing  correction  factor  of 

31 

1.32  for  the  impedance  of  the  effusion  orifice. 

Discussion 

The  principal  species  observed  mass  spectrometrically 

and  the  pressures  of  these  species  as  found  by  weight  loss 

rates  agree  well  with  the:  pressures  calculated  from  existing 

32 

thermodynamic  data.  Three  reactions  which  yield  vola¬ 
tile  species  in  excess  of  those  found  over  alumina  under 
neutral  conditions  are: 

a.  A1203  +  Nb  =  Al^  0(g)  +  Nb02(g) 

b.  A1203  +  3Nb  -  2 Al(g)  +  3Nb0(s ) 

c.  A1203  +  2Nb  +  1/2  02  =  Nb0(s)  +  Al20(g)  +  Nb02(g) 

Reaction  (c)  is  the  sum  of  reaction  (a)  and  the  oxidation 
of  niobium  to  form  niobium  oxide(s).  Reaction  (c)  is  limited 
by  the  oxygen  available  over  niobium  +  niobium  oxide(s)  as 
shown  in  curve  C  of  Figure  27.  Figure28  shows  the  calculated 
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partial  pressures  over  equations  (b)  and  (c);  equation  (a) 
is  equivalent  to  (c)  when  (c)  is  limited  by  oxvgen  over 
niobium  +  oxide  as  calculated  in  Figure  27.  The 
temperature  dependence  of  weight  loss  predicted  by  equa¬ 
tions  (b)  and  (c)  is  identical  in  slope  and  similar  in 
magnitude  to  the  observed  weight  loss  (Figure  26).  The 
predicted  weight  losses  are  twice  those  observed  at  all 
temperatures  studied. 

The  observed  weight  loss  rates  (Figure26)  are  about  two 

orders  of  magnitude  greater  than  evaporation  loss  of 

33 

alumina  under  neutral  conditions .  The  observed  reactions 

between  alumina  and  niobium  are  driven  by  the  stable 

volatile  products  Al(g)  or  ALO(g).  Calculations  using 

32,34 

available  thermodynamic  data  have  shown  that 
refractory  oxides  other  than  alumina  can  be  used  success¬ 
fully  with  niobium  because  their  gaseous  species  are  less 
stable.  Table  8  gives  calculated  gaseous  reaction  product 
pressures  over  some  likely  reactions  between  niobium  and 
some  oxides  at  2000°K.  Beryllia  weight  loss  will  be  about 
doubled  by  the  presence  of  niobium  and  niobium  oxide(s) 
will  be  formed  if  these  two  contact.  Yttria  (and  other  stable 
oxides  of  low  vapor  pressure  metals)will  not  react  appreciably 
with  niobium.  Table  8  shows  that  the  most  favored  yttria- 
niobium  reaction  yields  NbC^g)  and  YO(g);  the  total  pressure 
of  reaction  products  at  2000°K  is  about  an  order  of  magni¬ 
tude  less  than  the  oxide  decomposition  pressure.  This 
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TABLE  8. 


Predicted  Reactions  Between  Niobium  and  Some  Oxides  at  2000°K 


Oxide 

Oxide 

Decomposition 

Pressure 

(Atmosphere) 

Total  Pressure 
of  Reaction 

Reaction  Products 

A12°3 

0° 

i 

o 

A1203  +  3Nb  -  2Al(g)  +  3NbO(s) 
A1203  +  2Nb  +  1  /2  02  - 

io’6-2 

NbO(s)  +  Al20(g)  +  Nb02(g) 

l0-6. 6 

BeO 

io-7-4 

BeO  +  Nb  =  Be(g)  +  NbO(s) 

IO'7, 6 

Y2°3 

00 

00 

1 

o 

Y203  +  3Nb  =  2Y(r)  +  3NbO(s) 

1 o” 12 ■ 1 

Y203  +  Nb  2  Y0(g)  +  NbO(s) 

lO-iO.6 

2Y203  +  Nb  =  Nb02(g)  +  4  Y0(g ) 

IO*9,6 

Th02 

IO"9,  1 

Th02  +  Nb  ThO(g)  +  NbO(s) 

Th02  +  1/2  Nb  - 

io'16,2 

ThO(g)  +  1/2  Nb02(g) 

10’  1  °’  4 
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predicted  compatibility  and  the  similarity  between 
thermal  expansivities  of  niobium  and  them 

likely  candidate  materials  for  high  temperature  metal - 
ceramic  joins. 

Conclusions 

The  equilibrium  reaction  between  alumina  and  niobium 
preceeds  principally  as  Al^O^  +  3Nb  2Al(g)  +  3Nb0(s); 
a  secondary  reaction  of  almost  equal  importance  under 
neutral  conditions  is  A^O^  +  Nb  +  1 /2  -  Al^fHg)  + 

NbO^g)  +  NbO(s).  The  secondary  reaction  can  dominate 
under  oxidizing  conditions.  Weight  loss  due  to  the 
gaseous  reaction  products  of  the  two  reactions  is  about 
one  hundred  times  the  weight  loss  due  to  evaporation 
of  alumina  under  neutral  conditions.  and  some 

other  oxides  can  be  used  at  high  temperatures  with 
niobium  with  no  significant  reaction. 
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APPENDIX  A 


Method  for  Determining  the  Diffusivits  of  Oxygen  in  Mo  from  the  Average 
concentration  of  Rods 

1 

It  was  mentioned  previously  that  Fick's  second  law  diffusion  equation  can 

be  solved  for  an  infinitely  long  solid  cylinder  if  two  boundary  conditions  are 

2 

known  and  if  the  diffusivity  is  not  a  function  of  concentration.  The  solution 
of  the  diffusion  equation  for  the  above  conditions  is: 

'«  -Di  2t 

[ C(r )  -  C.]/[C  -  C.]  =  1-2/a  Z  e  n  J  (r,  -  )/J.(a,q  )  (A-l) 

'  '  o  iJ  .  o’  n'  1  ’  n' 

n=  1 

where  C(r)  is  the  concentration  of  oxygen  at  a  radius  r  after  time  t;  Ch  is 
the  initial  uniform  concentration  of  oxygen  in  the  cylinder;  C()  is  the  con¬ 
centration  of  oxygen  at  the  surface  of  the  rod;  a  is  the  rod  radius;  D  is  the 
diffusivity;  t  is  the  time;  r  is  the  radial  distance  from  the  center  of  the  rod; 
Jq  and  Jj  are  Bissel  functions  of  order  zero  and  one,  respectively;  and  an 

is  the  root  of  J  ,  a)  =  0.  If  the  solution  for  the  cylinder  is  written  in 
o  n  7 

2 

terms  of  two  dimensionless  parameters,  Dt/a  and  r/a,  a  family  of  curves 

is  obtained  when  [C{r)  -  C^]/[Cq  -  C.)  is  plotted  versus  r/a  (one  curve  for 

2 

each  value  of  Dt/a  ).  The  particular  family  of  curves  which  one  obtains 
depends  on  the  rate  of  supply  of  oxygen  to  the  external  cylinder  surface 
which  is  characterized  by  a  transfer  coefficient  3.  This  coefficient  is 

3 

analyogous  to  the  surface  coefficient,  HC,  applied  to  heat  flow  problems. 
The  remainder  of  the  discussion  in  this  appendix  applies  equally  well  to  a 
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family  of  curves  for  any  given  8.  No  two  curves  of  a  family  cross;  thus, 

2 

the  area  under  each  curve  is  uniquely  determined  by  Dt/a  .  The  area  under 

each  curve  of  the  family  of  curves  r  [C(r)  -  (T]/[Cq  -  C.]  vs  r/a  is  also 

2 

uniquely  determined  by  Dt/a  . 


The  area,  A,  under  each  curve  is  given  by, 


r  /a  =  1 

a  =  r 

A  J, 

r  /a  =  o 


(C(r)  -  C.) 

r  (C  -  C.) 
o  1 


d 


(A  -2 ) 


which  simplifies  to, 


A  = 


a(C  -C.) 
'  o  r 


r /a  =  1 

j*  C(r)  rdr 

r  /a  =  o 


r  /a  -  1 

C.  r"  rdr 

1  J  / 

r/a  =  o 


(A-3) 


2 

Multiplying  the  first  integral  of  equation  (A-3)  by  2TT't/2'T-l  (where  l  is  the 
cylindrical  rod  length)  and  integrating  the  second  integral,  we  obtain 


r  /a  =  1 

C(r)  £2nrdr 

r  /a  =  o 


C. 

l 


(A  -4) 


A  cylindrical  element  of  volume  of  thickness,  dr,  at  constant  radius,  r,  is 

given  by  2rr-trdr.  C(r)  is  the  oxyben  concentration  throughout  this  volume 

element;  thus  C(r)  2n-trdr  is  the  number  of  oxygen  atoms  in  the  volume  element, 
r  /a  -  o 

Thus,  J  C(r)  't2Trdr  equals  the  total  number  of  oxygen  atoms,  N,  in  the 

r  /a  =  1 

whole  cylinder  of  radius,  a,  L if  C(r)  has  units  of  oxygen  atomx/unit  volume]. 
Substituting  N  into  equation  (A -4)  we  find, 
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A 


1 


1 

2nl 


N  -  C. 
1 


(A-5) 


Consider  now  only  the  first  term  inside  the  brackets  of  equation  (A-5). 
Multiplying  the  top  and  bottom  of  this  term  by  the  total  volume,  V,  we  obtain 
V/2rr^’N/V;  but  N/V  is  the  average  concentration  of  oxygen,  C  ,  in  the 
rod.  Substituting  V/2n'L*C  into  equation  (A  -4)  and  remembering  V  =  na^-t 
we  obtain, 


A  = 


1 


2 

na  l 

2nl  av 


(A-6) 


which  simplifies  to 


A  = 


a(C  -  C.) 
av  i‘ 

2(C  -  C.) 

o  l' 


(A-7) 


Thus,  the  area  under  a  given  curve  corresponding  to  a  given  value  of  Dt/a 
is  related  to  the  average  concentration  of  oxygen  in  the  rod  by  equation 
(A-7)  and  C(r)  need  not  be  evaluated. 


The  most  useful  form  of  equation  (A-7)  for  the  present  work  is: 


C 

av 


2(C 


-Ci> 


C. 

l 


(A-8) 


92 


GEST -2062 


Appendix  A  Bibliography 

1.  "Technical  Summary  Report  of  Research  and  Development  Program 
of  Thermionic  Conversion  of  Heat  to  Electricity,"  Semi-Annual 
Report  GEST-2046, 

2.  J.  Crank,  "The  Mathematics  of  Diffusion,"  Oxford  at  the  Clarendon 
Press,  1956,  p.  66. 

3.  J.  B.  Austin,  "The  Flow  of  Heat  in  Metals,"  ASM  Book,  1942,  p.  120. 


GEST -2  062 


APPENDIX  B 

Estimation  of  Equilibrium  Solubility  of  Oxygen  in  Molybdenum  as  a  Function 
of  Temperature 

To  calculate  values  for  diffusivities  of  oxygon  in  molybdenum  from  Fick's 
first  law  at  temperatures  in  the  range  1500  to  1800°C,  the  oxygon  concentra¬ 
tion  gradients  in  the  molybdenum  had  to  be  known.  From  the  present  experi¬ 
mental  studies  with  molybdenum  capsules,  the  oxygon  concentration  on  the 
outside  capsule  surface  at  1800°C  was  found  to  be  5  ppm  and  the  concentra¬ 
tion  on  the  inside  capsule  surface  was  found  to  bo  39  ppm  at  1500°C.  Since 
this  data  was  insufficient  for  calculating  at  any  of  the  temperatures  of 
interest,  several  assumptions  were  made  to  obtain  approximate  values  of 
AC.  First,  the  outside  concentration  of  5  ppm  found  at  J800°C  was  used  for 
capsules  at  all  temperatures.  Second,  the  concentration  at  the  inside  surface 
of  the  capsule  at  each  temperature  was  estimated  from  the  experimental 

value  obtained  at  1500°C  and  oxygon  solubility  data  reported  by  Few  and 
1 

Manning.  They  reported  a  solubility  limit  for  oxygen  in  molybdenum  as  a 
function  of  temperature  as  shown  by  the  solid  portion  of  curve  A  of  Figure  B-l. 
The  dashed  portion  of  the  curve,  A,  is  an  extrapolation  of  their  data  to  1800°C. 
They  exposed  samples  to  oxygen  pressures  well  in  excess  of  equilibrium 
pressures  for  the  formation  of  molybdenum  oxygen  (solid)  and  determined 
solubility  limits  by  the  microscopic  metal  presence  or  absence  of  oxide 
in  their  samples.  This  procedure  and  analytical  technique  could  easily 
lead  to  high  values  of  the  solubility  limit.  The  oxygen  pressure  of  the 


94 


Temperature- 


GEST  -2062 


Curve  A  Few  and  Manning 

Curve  B  Interpolation  between  A  and  C 

Curve  C  Extrapolation  of  present  data 


FIGURE  B-l.  SOLUBILITY  LIMIT  OF  OXYGEN  IN  MOLYBDENUM 
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present  experiments  was  very  close  to  the  equilibrium  value.  Because  of 
this,  the  value  of  39  ppm  at  1500°C  is  believed  to  be  more  accurate  than 
Few  and  Manning's  data.  To  estimate  the  concentrations  at  the  inside 
capsule  surface  for  temperatures  higher  than  1500°C,  the  slope  of  Curve  A 
is  assumed  to  be  correct  and  a  curve  C  parallel  to  curve  A  and  passing 
through  39  ppm  at  1800°C  was  constructed.  In  addition  to  the  temperature 
dependence  shown  by  curve  C  we  would  expect  (from  Henry's  Law)  a  contri¬ 
bution  to  the  solubility  at  each  temperature  above  1500°C  from  the  increased 
equilibrium  pressure  relative  to  the  pressures  used  by  Few  and  Manning. 

The  pressures  used  by  Few  and  Manning  differed  from  the  equilibrium 
pressures  used  in  the  present  investigation,  as  shown  in  Table  B-I.  The 
pressure  ratio  (Table  B-I)  decreased  with  increasing  temperature  showing 
that  the  equilibrium  pressure  at  1800°C  was  approximately  one  order  of 
magnitude  or  25%  closer  to  the  pressure  used  by  Few  and  Manning  than  the 
equilibrium  pressure  at  1500°C.  The  pressure  contribution  to  the  solubility 
of  oxygen  at  1800°C  was  therefore  estimated  to  be  approximately  1/4  of  the 
concentration  difference  between  curves  A  and  C  or  point  (1)  of  Figure  B-l. 
The  solid  solubility  of  oxygen  at  intermediate  temperature  was  estimated 
by  constructing  the  dashed  straight  line  B  from  39  ppm  at  1500°C  to  53  ppm 
(point  (1))  at  1800°C.  The  resulting  solubility  limits  were  48,  44,  and  34 
ppm  at  1700,  1600,  and  1400°C,  respectively. 

3 

In  preliminary  experiments,  Rudman  found  a  solubility  of  oxygen  in  molyb¬ 
denum  of  45  f  5  ppm  at  1800°C.  This  tends  to  support  the  method  of  esti¬ 
mation  used  above  rather  than  the  data  of  Few  and  Manning. 
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TABLE  B-l 


’cst  Pressures 


Temperature 


1500 


Pressure  (microns) 

'1)  Few  &  Manning  (2)  Prosont  Stud1) 


2. 5  x  10' 


1 . 2  x  i  0- 


Pressure 

Ratio 

(1)  :  (2) 

2.  08  x  lO*1 


7.4  x  10' 


7.  5  x  I  O' 


i  x  10 


1. 9  x  10- 


3.  8  x  1  O' 


5xJ0 


4.  1  x  10- 


2. 2  x  10 


^Extrapolated  from  Few  and  Manning's  work. 
2 

**King  and  Weller  data  were  used. 
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